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ABSTRACT 

We present a catalog of revised effective temperatures for stars observed in long-cadence mode in the Kepler 
Input Catalog (KIC). We use SDSS griz filters tied to the fundamental temperature scale. Polynomials for 
griz color-temperature relations are presented, along with correction terms for surface gravity effects, metal- 
licity, and statistical corrections for binary companions or blending. We compare our temperature scale to the 
published infrared flux method (IRFM) scale for VtJKs in both open clusters and the Kepler fields. We find 
good agreement overall, with some deviations between (J - Ks)-hased temperatures from the IRFM and both 
SDSS filter and other diagnostic IRFM color-temperature relationships above 6000 K. For field dwarfs we find 
a mean shift towards hotter temperatures relative to the KIC, of order 215 K, in the regime where the IRFM 
scale is well-defined (4000 K to 6500 K). This change is of comparable magnitude in both color systems and in 
spectroscopy for stars with Teff below 6000 K. Systematic differences between temperature estimators appear 
for hotter stars, and we define corrections to put the SDSS temperatures on the IRFM scale for them. When 
the theoretical dependence on gravity is accounted for we find a similar temperature scale offset between the 
fundamental and KIC scales for giants. We demonstrate that statistical corrections to color-based tempera- 
tures from binaries are significant. Typical errors, mostly from uncertainties in extinction, are of order 100 K. 
Implications for other applications of the KIC are discussed. 

Subject headings: stars: fundamental parameters 



1. INTRODUCTION 

One of the most powerful applications of stellar multi-color 
photometry is the ability to precisely infer crucial global prop- 
erties. Photometric techniques are especially efficient for 
characterizing large samples and providing basic constraints 
for more detailed spectroscopic studies. Modern surveys fre- 
quently usedfilters designed for the Sloan Digital Sky Survey 
(SDSS: lAihara e t al. 2011), however, while traditional corre- 
lations between color and effective temperature (Teff), metal- 
licity ([Fe/H]), and surface gravity (logg) have employed 
ot her filter sets, ty pically on the Johnson-Cousins system. 
In lAn et al.1 (l2009ai hereafter A09) w e used SP SS photom- 
etry of a solar-metalHcity cluster M67 (lAn et al.l [2b08) to de- 
fine a photometric ugriz-T^n relation, and checked the metal- 
licity scale using star clusters over a wide range of metal- 
licity. This scale wa s applied to the Virgo overdensity in 
the halo by I An et al.l (l2009bi) . The approach used is sim- 
ilar in s pirit to earher wo rk in the Johnson-Cousins filter 
system dPinsonneaulietiD l200l l2004t lAn et all l2007allbh : 
the latter ef fort use d the c olor-temperature relationships of 
iLeieune et al.l 019971 119981) with empirical corrections based 
on cluster studies. 

A revised color-temperature-metallicity relationship 
for late-type stars has recently been published by 
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ICasagrande et al.l (120101 hereafter CIO); it is based on 
the infrared flux method (IRFM). There are a number of 
advantages of this approach, as discussed in CIO, but there 
is a lack of native SDSS data in the stars used to define 
the calibration itself. Fortunately, the color- temperature 
relationships in CIO are defined for JH Ks colors in the Two 
Micron All Sky Survey (2MASS; Ski-ut skie et al.l l2006). and 
the Kepler mission provides a large body of hig h quality griz 
photo metry for stars in the 2MASS catalog (iBrown et alJ 
l20Tlh . 

In this paper we use griz data in the Kepler Input Cata- 
log (KIC) in conjunction with 2MASS to compare the effec- 
tive temperature scale for the griz colors to the IRFM scale. 
For this initial paper we concentrate on the mean relation- 
ships between the two systems for the average metallicity of 
the field sample, taking advantage of the weak metallicity 
dependence of the color- Teff relationships that we have cho- 
sen. In a follow-up paper we add information from spectro- 
scopic metallicity and \ogg determinations to compare empir- 
ical photometric relationships involving these quantities to the 
theoretical relationships used in the current work. Unresolved 
binaries and extinction errors can be severe problems for pho- 
tometric temperature estimates, and another goal of this work 
is to quantify their importance. 

Another important matter, which we uncovered in the 
course of our research, concerns systematic errors in the griz 
photometry in the KIC. For large photometric data sets, it can 
be difficult to assess such errors. Fortunately, we can also 
compare photometry used in the KIC with photometry in the 
same fields from the SDSS; the latter is important for numer- 
ous applications of data derived from the Kepler mission. We 
will demonstrate that there are significant systematic differ- 
ences between the two, and derive corrections to minimize 
these effects. 

We therefore begin with a discussion of our method in Sec- 
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tion|2] Alon g with a description on the sample selection in the 
KIC (§ 12. Il l, we compare the SDSS and KIC ph otom etry and 
derive corrected KIC magnitudes and colors (§ I2.2l i. A ba - 
sis model isochrone in the SDSS colors is presented (§ l2.3l l, 
and a met hod o f determining photometric T^ff from griz is de- 
scribed (§ iMli. Both the IRFM/VJHK, and SDSS/griz tem- 
perature scales are compared to the KIC dwarf temperatures 
in Section[3] where a ~ 200 K offset is found in the KIC with 
respect to both IRFM and SDSS temperature scales. We also 
present a method of correcting the dwarf temperature scale for 
giants (§ 13. 2t . For well-studied open clusters, we find a good 
agreement overall between SDSS and IRFM, but find some 
systematic deviations between IRFM (J - Ks)-hased tempera- 
tures from the IRFM and both SDSS filte r an d other diagnos- 
tic IRFM color-temperature relations (§ I3.3l l. We provide a 
formula to put SDSS Teff on the consistent scale with IRFM. 
These findings are confirmed using spectroscopic temperature 
determinations (§ |3.4t . We also discuss the impact of unre- 
solved binari es a nd uncertainties in the extinction estimates 
(§ 13.51 and § 13.6b . Our revised catalog is presented in Sec- 
tion H] where we provide a recipe for estimating Teff for inter- 
ested readers, if the application of our technique is desired to 
the entire KIC sample in general. We discuss the implications 
of our new fundamental Teff scale in Section|5] 

2. METHOD 

Our basic data come from the long-cadence sample in the 
KIC. From this we extracted a primary sample of dwarfs in 
the temperature range where our calibrations are best con- 
strained; our procedure is given in Section lZTI We uncovered 
some offsets between KIC and nati ve SD SS photometry, and 
describe correction terms in Section l2!2l Our methods for de- 
riving co lor-t em pera ture relationships in griz are described in 
Sections |23] and IZ4l 

2.1. Sample 

We took griz photometry from the KIC (iBrown et al.l201 ll) : 
photometric uncertainties were taken as 0.01 mag in gri and 
0.03 mag in z- Errors were taken from the quadrature sum of 
uncertainties in the individual filters. JHKs photometry was 
taken from the All Sk y Data Release of the 2MASS Point 
Source Catalog (PSC: ISkrutskie et all l2006fL and checked 
against complementary information in the KIC itself. 

For our sample we chose long-cadence targets in the KIC; 
our initial source had 161,994 candidates. We selected stars 
with griz photometry detected in all of the bandpasses. This 
sample was nearly complete in the 2MASS catalog. We ex- 
cluded a small number of sources with 2MASS photometry 
quality flags not equal to AAA (A^ = 3,602) and stars with 
colors outside the range of validity of either the IRFM or 
SDSS scales (A^ = 11,830), leaving us with a main sample 
of 146,562 stars. We then further restricted our sample by 
excluding stars with log^ estimates below 3.5 dex in the KIC 
(A^= 19,663) for a dwarf comparison sample of 126,899. We 
illustrate the distribution of stars in the sample in 100 K bins 
in Figure [U both in the initial catalog (top panel) and the re- 
vised one in this paper (bottom panel). We did not use the gi- 
ants in our comparison of the dwarf -based temperature scale 
(Section |3]l, but we do employ theoretical \ogg corrections 
to the photometric temperatures for the purposes of the main 
catalog (see Sections [3. 2l and l4r2] i. 

^ See |http://www.ipac.caltech.edu/2mass7] 
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Fig. 1. — Long-cadence data from the KIC (top) and with our revised 
SDSS-based effective temperature scale (bottom). Data are binned in 100 K 
increment. Dwarfs with KIC logg > 3.5 (open histogram) are separated from 
giants with lower logg (shaded histogram). 



2.2. Recalibration of the KIC Photometry 

We adopted three primary color indices (g - r, g - i, and 
g-z)as our temperature indicators for the SDSS filter system. 
A preliminary comparison of colors yielded surprising inter- 
nal differences and trends as a function of mean Teff in the 
relative temperatures inferred from these color indices (see 
below). Because the A09 color-color trends were calibrated 
using SDSS photometry of M67, this reflects a zero-point dif- 
ference between the KIC and SDSS photometry in the color- 
color plane. It is not likely that this difference is caused by 
extinction or stellar population differences because all three 
colors have similar sensitivities to extinction and metallicity. 
Initi ally, we suspec ted problems with the SDSS calibration 
(see lAn et al.ll2008l for a discussion of zero-point uncertain- 
ties). However, the differences seen were outside of the error 
bounds for the SDSS photometry. For a fraction of the tar- 
gets (about 2%) the temperatures inferred from different color 
sources (SDSS versus IRFM from 2MASS colors) are also 
discordant by more than three standard deviations, in some 
cases by thousands of degrees in Tgff. We examine both phe- 
nomena below. 

About 10% of the stars in the Kepler field are covered in the 
most recent data re lease (DR8) of the SDSS imaging survey 
jAihara et alj|201 ll) . There is an overlap in the two photomet- 
ric sets at 14<r< 18. We compare photometry for stars in 
common in Figure |2l With a 1" search radius, we found that 
the median differences (in the sense of the SDSS minus KIC), 
after rejecting stars with differences greater than 0.2 mag on 
both sides, are Ag = -0.040, Ar = -0.028, A/ = -0.045, and 
Az = -0.042. 

Inspection of Figure|2]shows that these differences are also 
functions of color. Solid lines are a linear fit to the data after 
an iterative 3a rejection. The linear transformation equations 
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Fig. 2. — Photometry comparisons between SDSS (DR8) and the KIC in 
the sense of the former minus the latter. Comparisons are shown in griz from 
top to bottom panels. Solid lines are a linear fit to the residuals. 

are as follows: 

^SDSS = ^Kic + 0.092 1 (g - r)Kic - 0.0985 , ( 1 ) 

rsDSS = rKic + 0.0548(r- Okic" 0.0383, (2) 

iSDSS = !Kic + 0.0696(r- Okic -0.0583, (3) 

zsDSS =zkic + 0.1587(/-z)kic-0.0597, (4) 

where the subscripts indicate either SDSS or KIC photometry. 

It is possible that the SDSS photometry in the Kepler field 
has some zero-point shifts with respect to the main SDSS sur- 
vey database; the SDSS photometry pipeline can fail to work 
properly if the source density is too high. To check this, we 
compared the KIC and the DAOPHOT crowded-field photom- 
etry (.Anetal...2008.) in the NGC 6791 field. Although the 



sample size is smaller, a comparison of DAOPHOT and KIC 
photometry in the NGC 6791 field yields systematic offsets 
in the same sense as the field mean in all band passes. Given 
that the cluster fiducial sequence from DAOPHOT photome- 
try matches that from an independent study ( Clem et al.i2008b 
relatively well (A n et a l. 2008), it is unlikely that the offsets 
seen in Figure|2]are due to zero-point issues in the SDSS pho- 
tometry. 

We also checked the standard star photometry in 
iBrown et alj (12011'), which is originally from the SDSS DRl 
photometry for 284 stars outside of the Kepler field. We com- 
pared with the SDSS DR8 photometry, but did not find the 
aforementioned trends outside of those expected from ran- 
dom photometric errors: the mean differences (SDSS minus 
KIC values) were -0.009, -0.004, +0.004, and +0.012 mag 
in griz, respectively, with an error in the mean of the order 
of 0.001 mag. Therefore, revisions of the standard magni- 
tude system (SDSS DRl versus DR8) do not appear to be the 
explanation either We also investigated the possibility of a 
zero-point difference between the faint and bright stars in the 
KIC, which had different exposures. However, we found that 
magnitude offsets with respect to SDSS are similar for both 
samples, and that the internal dispersion of the KIC tempera- 
ture estimates is essentially the same with each other Regard- 
less of the origin, the differences between the SDSS and KIC 
photometry are present in the overlap sample, and we there- 
fore adjusted the mean photometry to be on the most recent 
SDSS scale. 

Inspection of Figure |2] also reveals another problem in the 
KIC photometry: a sub-population of stars are much brighter 
in the KIC than in the SDSS even after photometric zero-point 
shifts have been accounted for. We do believe that unresolved 
background stars explains the occasional cases where differ- 
ent colors predict very different temperatures. In Figure |2] 
there are many data points that have KIC magnitudes brighter 
than the SDSS ones. We attribute these stars to blended 
sources in the KIC. The mean FWHM of SDSS images is 
1 .4", while that of KIC photometry is 2.4". 

To check on this possibility, we cross-checked 13,284 stars 
in common between DR8 of the SDSS and our KIC sample. 
312 stars had a resolved SDSS source within 2.4", while 20 
have two or more such blended sources. 2.5% of the stars 
would therefore have resolved blends between the resolution 
of the two surveys. If we assume that the space density of 
blends is constant, we can use the density of blends to estimate 
the fraction present even in the higher resolution SDSS sam- 
ple. When this effect is accounted for, we would expect 3.8% 
of the KIC sources to have a blended star within the resolution 
limit of the KIC. The average such star was 2.85 mag fainter 
than the KIC target, sufficient to cause a significant anomaly 
in the inferred color- temperature relationships. A comparable 
fraction of the catalog is likely to have similar issues. A sig- 
nificant contribution from background stars would in general 
combine light from stars with different temperatures. As a 
result, one would expect different color-temperature relations 
to predict discordant values. We therefore assess the internal 
consistency of the photometric temperatures as a quality con- 
trol che ck in our revised catalog to identify possible blends 
(Sectionl24b. 

To identify blended sources in the KIC, we further per - 
formed a test using the separation index dStetson et alj|2003l) . 
which is defined as the logarithmic ratio of the surface bright- 
ness of a star to the summed brightness from all neighboring 
stars (see also .An et al.„2008i) . However, we found that apply- 
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ing the separation index to the KIC does not necessarily pro- 
vide unique information for assessing the effects of the source 
blending. 

2.3. Base Model Isochrone 

We adopted stellar isochrones in A09 for the estimation of 
photometric temperatures. Interior models were computed 
using YREC, and theoretical color-Teff relations were de- 
rive d from th e MARC S stellar atmospheres model: see A09 
and lAn et al.l (l2009bh for details. These model colors were 
then calibrated using observed M67 sequences as in our ear- 
lier work in the John son-Cousins system ( Pi nsonneault et^ 
I2003L I2004t lAn et al.. 2007a b). The empirical color correc- 
tions in ugriz were defined using M67 at its solar metallicity, 
and a linear ramp in [Fe/H] was adopted so that the color cor- 
rections become zero at or below [Fe/H]< -0.8. Detailed test 
on the empirical color corrections will be presented elsewhere 
(An et al. 2012, in preparation). 

As a base case of this work, we adopted the mean metal- 
licity recorded in the KJC of [Fe/H]= -0.2. This metal- 
licity is comparable to, or slightly below, that in the so- 
lar neighborho od. For e x ample , the Geneva-Copenhagen 
Survey (Nordstrom et al. I 120041) has a mean [Fe/H] of 
—0.14 dex with a disper sion of 0.19 dex; a recent revision by 
ICasagrande et all (1201 ll) raises the mean [Fe/H] to —0.07 dex, 
which is a fair reflection of the systematic uncertainties. The 
bulk of the KIC dwarfs are about 100 pc above the galac- 
tic plane, and thus would be expected to have somewhat 
lower metallicity. In the following analysis, we assumed 
[Fe/H]= -0.2 when using griz- or IRFM color- Teff relation- 
ships, unless otherwise stated. 

Table [T] shows our base model isochrone at [Fe/H]= -0.2 
and the age of 1 Gyr All colors are color-calibrated as de- 
scribed above. Note that the isochrone calibration is defined 
for the main-sequence only; the relevant corrections for the 
lower gravities of evolved stars are described separately in 
Section [32] The SDSS photometry did not cover the main- 
sequence turn-off region of M67 because of the brightness 
limit in the SDSS imaging survey at r ~ 14 mag. As a re- 
sult, the M67-based griz color calibration is strictly valid at 
4000 < Teff < 6000 K (see Figure 17 in A09). 

The choice of 1 Gyr age in our base model isochrone has 
a negligible effect on the color- Teff relations. The difference 
between 1 Gyr and 12 Gyr isochrones is only less than 5 K 
near main-sequence turn-off. However, younger age of the 
models enables the determination of photometric Teff over a 
wider range of colors at the hot Teff end. 

From Table [T] we derived polynomial color- Teff relations 
of our base model for convenience of use. The following 
relationship was used over the temperature range 4080 K 
<reff(YREC)<7000/:: 

5040/reff = aQ + aix+ a2X^ + a^x^ + a4X^ + asx^ (5) 

where x represents g - r, g - i, or g - z, and ao-a^ are coef- 
ficients for each color index as listed in Table |2] Difference 
in Teff inferred from these polynomial equations compared to 
those found in Table [T]from interpolation in the full tables are 
at or below the 6 K level. 

In Table |3] we provide the metallicity sensitivity of the 
color- Teff relations in the model isochrones at several [Fe/H]. 
To generate this table, we compared 1 Gyr old isochrones 
at individual [Fe/H] with our fiducial model (Table [U at 
[Fe/H]= -0.2 for each color index, and estimated the Teft dif- 
ference at a given color (individual models minus the fiducial 



isochrone). We include the sensitivity to metallicity predicted 
by atmospheres models, but do not include an additional em- 
pirical correction below [Fe/H]= -0.8 because the cluster data 
did not require one. The Teff at a fixed color generally be- 
comes cooler at a lower [Fe/H]. We use the metallicity cor- 
rections in the comparisons with spectroscopic T^ff w here we 
have reliable [Fe/H] measurements (see Section I3.4l l , but do 
not apply corrections to the KIC sample (see Sections [3. H and 

2.4. Photometric r^ff Estimation 

The stellar paramete rs for the KIC wer e generated using a 
Bayesian method (see iBrown et alj[201 1 , for a discussion). 
We adopt a less ambitious approach focused on KIC stars 
identified as dwarfs. The three key assumptions in our work 
are that we define Tgff at a reference [Fe/H] and the model 
logg (Table [T]), and that we adopt the map-based E{B - V) in 
the KIC as a prior. Within this framework we can then derive 
independent temperature estimates from the griz photometry 
and infer the random Teff errors. Uncertainties in the extinc- 
tion, the impact on the colors of unresolved binaries, and pop- 
ulation (metallicity and logg) differences can then be treated 
as error sources. In the latter case, we can compute correction 
terms to be used if there is an independent method of mea- 
surement. This approach is not the same as the one that we 
have employed in earlier studies, so a brief justification is in 
order. 

The traditional approach to photometric parameter estima- 
tion is to take advantage of the fact that different filter com- 
binations respond to changes in metallicity and extinction. If 
one has the proper template metallicity and extinction, for ex- 
ample, the answers from the various colors will agree within 
photometric errors; if not, the pattern of d ifferences can be 
used to solve for them (see lAn et al.ll2007"allbh . 

The particular problem for the KIC is that the available 
color combinations in griz are rather insensitive to both over 
the narrow metallicity range and the m odest mea n extinctions 
(0 < E{B -V)< 0.2) in the field (see lAn et al.l l2b09b, for a 
discussion of gr/z-based estimates). In other words, all the 
color combinations in griz produce similar metallicity sen- 
sitivities of color- Teff relations. Therefore, even though the 
absolute change of photometric Teff can be significant by the 
error in the adopted metallicity, it is difficult to infer photo- 
metric metallicities based on the available filter combinations 
in griz alone. 

The temperature estimates in iLejeune et alJ (119971 119981) . 
which were used as the basis color- Teff relations in our prior 
color calibration in the Johnson-Cousins system (lAn et alJ 
l2007allb ). are insensitive to logg near the main sequence, and 
the IRFM scale in CIO does not include an explicit logg de- 
pendence for the temperatures. As a result, we believe that 
the most fruitful approach is to define a benchmark tempera- 
ture estimate. If additional color information or spectroscopic 
[Fe/H] data become available, the relevant corrections can be 
applied, and we present methods below to do so (Section[3]). 

The KIC gravities for cool stars are precise enough to sep- 
arate dwarfs (KIC logg > 3.5) from giants (KIC logg < 3.5) 
and to be used as a basi s for corrections to the temperatures 
for giant stars (Sections 13.21 and l 4.2b . The KIC metallicities 
are more problematic, and we do not use them for temperature 
corrections. Instead the metallicity sensitivity is included as 
an error source in our effective temperature estimates. 

We adopted the map-based KIC catalog extinction esti- 
mates (Ay) and the .Cardelli et al.. (1989) standard extinction 
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0.593 


0.785 


0.869 


0.82 


5207.7 


-0.44 


4.61 


5.80 


0.619 


0.822 


0.914 


0.81 


5142.1 


-0.48 


4.62 


5.90 


0.647 


0.861 


0.961 


0.79 


5077.2 


-0.51 


4.63 


6.00 


0.676 


0.901 


1.010 


0.78 


5013.1 


-0.55 


4.63 


6.10 


0.707 


0.943 


1.061 


0.77 


4949.9 


-0.58 


4.64 


6.20 


0.739 


0.987 


1.115 


0.75 


4887.6 


-0.62 


4.64 


6.30 


0.773 


1.034 


1.170 


0.74 


4826.3 


-0.65 


4.65 


6.40 


0.808 


1.082 


1.228 


0.73 


4766.1 


-0.68 


4.65 


6.50 


0.846 


1.132 


1.288 


0.72 


4707.4 


-0.72 


4.66 


6.60 


0.885 


1.185 


1.350 


0.71 


4650.3 


-0.75 


4.66 


6.70 


0.926 


1.239 


1.414 


0.70 


4595.1 


-0.78 


4.66 


6.80 


0.967 


1.295 


1.479 


0.69 


4541.8 


-0.81 


4.67 


6.90 


1.007 


1.350 


1.544 


0.68 


4490.6 


-0.84 


4.67 


7.00 


1.047 


1.406 


1.610 


0.67 


4441.5 


-0.87 


4.68 


7.10 


1.084 


1.461 


1.675 


0.66 


4394.5 


-0.90 


4.68 


7.20 


1.121 


1.515 


1.740 


0.65 


4349.3 


-0.93 


4.69 


7.30 


1.155 


1.568 


1.804 


0.64 


4306.0 


-0.96 


4.69 


7.40 


1.188 


1.619 


1.867 


0.63 


4264.4 


-0.99 


4.70 


7.50 


1.218 


1.670 


1.928 


0.62 


4224.3 


-1.02 


4.70 


7.60 


1.247 


1.719 


1.989 


0.61 


4185.9 


-1.04 


4.71 


7.70 


1.273 


1.766 


2.048 


0.60 


4148.9 


-1.07 


4.71 


7.80 


1.298 


1.812 


2.105 


0.59 


4113.7 


-1.10 


4.72 


7.90 


1.320 


1.856 


2.160 


0.58 


4079.8 


-1.13 


4.73 


8.00 


1.340 


1.898 


2.213 


0.58 


4047.1 


-1.16 


4.74 


8.10 


1.358 


1.939 


2.265 



curves with Ay = 3.1E(B — V). Extinction coefficients in griz 
were derived in A09: Ag = 1.196Ai/ , A^ = 0.814Av , Ai = 
0.672Ay, and A- = OASSAy. We further took Aj = 0.282Ay, 
Ah = O.lSOAy, Ak^ = OAllAy, an d Aw = 1 .OSOA y, where Vj 
represents the Tycho V passband (lAn et al.ll2007ai) . 

For a given extinction-corrected set of griz magnitudes, 
we searched the best-fitting stellar template in the model 
isochrone for each star in the KIC. The mean T^ff was ob- 
tained by simultaneously fitting the models in griz, assuming 
0.01 mag error in gri and 0.03 mag error in z- We also esti- 
mated Teff using the same model isochrone, but based on data 
from each of our fundamental color indices (g - r, g - i, and 



g — z), which is simply a photometric Teff estimation from a 
single color- Teff relation. Its purpose is to readily identify and 
quantify the internal consistency of our primary temperature 
determination from the multi-color- Teff space. 

In Figure |3] we plot the internal dispersion and the mean 
trends of Tgff from a given color index with respect to the aver- 
age error-weighted temperature from griz for all of the dwarfs 
in our sample. The top panel shows the case of the original 
KIC data, and the bottom panel shows the one for the cor- 
rected KI C ph otometry. The magnitude corrections described 
in Section l2!2l were motivated by concordance between SDSS 
and the KIC. Nevertheless, the results when using the recal- 
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TABLE 2 

Coefficients for Polynomial CoLOR-Teff Relations 



Coeff. 








fo 


0.6676 


0.6888 


0.7053 


ai 


0.3434 


0.2012 


0.2022 


02 


0.5851 


0.4518 


0.2733 


03 


-0.6919 


-0.4871 


-0.2844 


£24 


0.1445 


0.1926 


0.1079 


05 


0.0594 


-0.0256 


-0.0144 



Note. — Coefficients in equation [5] Tliese cofficients are valid at 4080 < T^-fr < 
7000 K, orO.I3 < - r)o < 1.34, 0.13 < - Oo < 1.90, and 0.07 < (g - z)o < 2.21. 
respectively. 



T 1 1 1 1 1 1 1 1 1 r 




I I I I I I I I I I I 

0.0 0.5 1.0 

g - r 

Fig. 3. — Internal dispersion in T^u estimates (solid line) and differences 
between the mean gnz-based T^ff and that inferred from g — r (long-dashed), 
g - i (short-dashed), and g — z (dotted). Original KIC photometry is used in 
the top panel, while corrected KIC photometry is used in the bottom panel. 

ibrated KIC photometry as temperature indicators were ex- 
tremely encouraging. 

Although the internal agreement is not complete, the re- 
maining differences in the bottom panel of Figure [3] are com- 
parable to the zero-point uncertainties discussed in lAn et alj 
(12008) . We view this as strong supporting evidence for the 
physical reality of the magnitude corrections illustrated in 
Figure |2] We therefore recommend that the zero-points of 
the KIC photometry be modified according to equations 1-4. 
In the remainder of the paper, we use magnitudes and colors 
adjusted using these equations. 

3. REVISED Teff SCALE FOR THE KIC 

We begin by evaluating the Tgff inferred from the IRFM and 
the SDSS systems for dwarfs (KIC \ogg > 3.5). We then 
use open clusters and comparisons with high-resolution spec- 
troscopy to establish agreement between the two scales, in- 
dicating the need for correction to the KIC effective tempera- 
tures. We then evaluate the impact of binaries, surface gravity. 



and metallicity on the colors. We provide statistical correc- 
tions to the temperatures caused by unresolved binary com- 
panions, as well as corrections for log^ and metallicity. We 
then perform a global error analysis including extinction un- 
certainties and the mild metallicity dependence of our color- 
temperature relationships. The latter is treated as a tempera- 
ture error source because we evaluate all KIC stars at a mean 
reference metallicity ([Fe/H]= -0.2). 

3.1. Temperature scale comparisons for dwarfs 

We have three native temperature scales to compare: the 
one in the KIC, our isochrone-based scale from griz (here- 
after SDSS or gr/z-based scale unless otherwise stated), and 
one from the (7 - A'i)-based IRFM. Below we compare the 
mean differences between them and compare the dispersions 
to those expected from random error sources alone. We find 
an offset between the KIC and the other two scales. The IRFM 
and SDSS scales are closer, but some systematic differences 
between them are also identified. In this section, we examine 
various effects that could be responsible for these differences, 
and finish with an overall evaluation of the error budget. 

We computed IRFM and SDSS T^if estimates assuming 
[Fe/H]= -0.2. In terms of the temperature zero-point, adopt- 
ing this metallicity led to mean shifts of +20 K in / - K^, 
and -40 K in the gr/z-based Teff estimate, relative to those 
which would have been obtained with solar abundance. In 
other words, changes in the adopted mean metallicity would 
cause zero-point shifts of 60 K in the overall Tgff scale com- 
parison. On the other hand, a scatter around the mean metal- 
licity in the Kepler field is another source of error that would 
make the observed Teff comparison broader. We discuss this 
in Sec tion [J!6l along with other sources of uncertainties. 

In the comparisons below we repeatedly clipped the sam- 
ples, rejecting stars with temperature estimates more than 
three standard deviations from the mean, until we achieved 
convergence. This typically involved excluding about 1 % 
of the sample. Such stars represent cases where the extinc- 
tion corrections break down or where the relative colors differ 
drastically from those expected for single unblended stars. 

Random errors were taken from the photometric errors 
alone and yield a minimum error in temperature. For the 
SDSS colors we also computed the internal dispersion in the 
three temperature estimates from g - r, g - i, and g - z, and 
used the larger of either this dispersion or the one induced 
by photometric errors as a random uncertainty. Median ran- 
dom errors for the SDSS and IRFM temperatures were 40 K 
and 171 K, respectively. These estimates are consistent with 
expectations from the observed dispersions of the colors (see 
Figure [TtTi. We then compared stars at fixed KIC temperature 
and computed the average Teff difference between those in- 
ferred from the IRFM, those inferred from griz, and the scale 
in the KIC itself. For a limited subset of stars, we also had 
Tyclio photometry, and computed temperatures from Vt - Ks. 
This sample is small, so we used it as a secondary temperature 
diagnostic. 

In Figures |4]-(5] we illustrate the differences between KIC 
and the IRFM and SDSS, respectively. For the IRFM scale in 
Figure m we compare Teff from J - K,. In Figure |5] we com- 
pare the mean SDSS temperatures inferred from griz to that 
in the KIC. In both cases we see a significant zero-point shift, 
indicating a discrepancy between the fundamental effective 
temperature scale and that adopted by the KIC. 

The Teff from the IRFM and the SDSS from individual color 
indices (g - r, g - i, and g - z) are compared in Figure |6] The 
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TABLE 3 

Teff Corrections for Different [Fe/H] 



[Fe/H] 



color 


-2.0 


-1.5 


-1.0 


-0.6 




-0.4 


-0.2" 


+0.0 


+0.2 


+0.4 


g-r 


0.1 


44 


36 


20 


19 




7 












0.2 


-32 


-32 


-27 


-6 




-4 





20 


61 


100 


0.3 


-120 


-103 


-74 


-33 




-20 





33 


84 


145 


0.4 


-162 


-144 


-110 


-60 




-32 





44 


104 


176 


0.5 


-168 


-165 


-130 


-75 




-42 





44 


116 


194 


0.6 


-186 


-185 


-144 


-80 




-42 





46 


113 


192 


0.7 


-214 


-195 


-149 


-84 




-44 





48 


104 


178 


0.8 


-238 


-198 


-149 


-85 




-42 





45 


97 


157 


0.9 


-262 


-204 


-153 


-86 




-43 





44 


83 


131 


1.0 


-282 


-212 


-160 


-94 




-46 





45 


77 


113 


1.1 


-289 


-223 


-174 


-107 




-51 





45 


69 


93 


1.2 




-233 


-207 


-136 




-60 





42 


45 


28 


1.3 












-83 





37 


-33 






0.1 


86 


67 


33 


20 




9 





-5 


-23 




0.3 


5 


-3 


-13 


-5 




-3 





16 


35 


61 


0.5 


-69 


-70 


-61 


-37 




-21 





27 


63 


108 


0.7 


-84 


-99 


-88 


-52 




-27 





33 


75 


127 


0.9 


-126 


-142 


-118 


-65 




-34 





35 


75 


125 


1.1 


-168 


-156 


-128 


-72 




-36 





39 


71 


112 


1.3 


-204 


-166 


-134 


-78 




^0 





41 


68 


99 


1.5 


-233 


-180 


-143 


-83 




-42 





44 


68 


95 


1.7 




-206 


-172 


-95 




-47 





50 


70 


93 


1.9 












-56 





61 


80 


101 


g-z 


0.1 


94 


68 


33 


15 




6 





-4 


-16 




0.3 


45 


25 


2 


-1 










10 


19 


33 


0.5 


-11 


-23 


-35 


-24 




-13 





20 


41 


72 


0.7 


-27 


-43 


-47 


-29 




-17 





20 


50 


87 


0.9 


-63 


-84 


-78 


-46 




-23 





26 


54 


92 


1.1 


-109 


-119 


-102 


-59 




-29 





31 


60 


93 


1.3 


-145 


-137 


-116 


-66 




-34 





35 


60 


90 


1.5 


-177 


-147 


-121 


-70 




-36 





39 


61 


87 


1.7 


-205 


-160 


-128 


-73 




-36 





44 


66 


92 


1.9 


-225 


-182 


-144 


-77 




-39 





49 


72 


97 


2.1 








-87 




-44 





58 


82 


109 



Note. — The .sense of the difference is the model Toff at a given [Fe/H] minus that of the fiducial metallicity, [Fe/H]= -0.2. The Tefr at a tixed color generally becomes cooler at a 
lower [Fe/H]. In other words, the above correction factor should be added to the SDSS 7[.fr, if the metallicity effects should be taken into account. 
" Fiducial metallicity. 



500 



-500 




4000 



5000 



6000 
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(KIC) 



7000 



in 

CO 
Q 
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Fig. 4. — Comparisons of the temperatures inferred from IRFIVI (J — Kg) 
as a function of KIC T^g. The color coding indicates the logarithmic number 
density of stars with a temperature and temperature difference at the indicated 
point (see legend). 



Fig. 5. — Same as in Figure|4] but from griz colors. 

IRFM scale for the Tycho Vj and 2MASS is used in the 
bottom right panel. The central result (that the KIC scale is 
too cool) is robust, and can also be seen in comparisons with 
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4000 5000 6000 7000 4000 5000 6000 7000 

Teff (KIC) (SDSS) 



Fig. 6. — Differences in Teff between tlie IRFM and SDSS scales as a function of KIC T^ff. T^ffig — r) vs. T^tdJ - Ks) (upper left); T^fiig — i) vs. T^niJ — Ks) 
(upper right); T^ffig — z) vs. T^ffiJ — Ks) (lower left); T^ffigriz) vs. rcff(^T — Ks) (lower right). The color coding defines the logarithmic number density of points 
with the indicated temperature and temperature difference (see legend for details). 



hig h-res olution spectrosco pic t emperature estimates (see Sec- 
tion |3.4| below). In Section |4~2l we provide quantitative tabular 
information on the statistical properties of the sample. 

The two fundamental scales (IRFM and SDSS) are close, 
but not identical, for cooler stars; they deviate from one an- 
other and the KIC above 6000 K (on the SDSS scale). As dis- 
cussed in Section [33] below. the total internal dispersion in the 
griz temperature estimates is also consistently larger for cool 
stars than that expected from random photometric uncertain- 
ties alone, and there are modest but real offsets between the 
two fundamental scales even for cool stars. We therefore need 
to understand the origin of these differences and to quantify 
the random and systematic uncertainties in our temperature 
estimates. 

Open clusters provide a good controlled environment for 
testing the concordance of the SDSS and IRFM scales. The 
SDSS scale was developed to be consistent with Johnson- 
Cousins-based temperature calibrations in open clusters, so 
a comparison of the An et al. and IRFM Johnson-Cousins 
systems in clusters will permit us to verify their underlying 
agreement. As we show below, the two scales are close for 
cool stars when B - V ,V - Ic, or V - Ks indices are employed 
in the temperature determinations, but exhibit modest but real 
systematics for the hotter stars. The IRFM relation in 7 - Ks, 
on the other hand, is found to have a systematic difference 
from those of these optical-2MASS indices. For the reasons 
discussed in the following section, we therefore adopt a cor- 



rection to our SDSS temperatures for hot stars, making the 
two photometric systems consistent. 

We can also check our methodology against spectroscopic 
temperature estimates, and need to consider uncertainties 
from extinction, binary companions, and metallicity. We 
therefore begin by defining an extension of our method to gi- 
ants, which can be checked against spectroscopy. We then 
look at open cluster tests, spectroscopic tests, binary effects, 
and the overall error budget. 

3.2. Tests of the temperature scale for giants 

Our YREC Teff estimates are based on calibrated isochrones 
(Table [Til, which do not include evolved stars. About 14% of 
the KIC sample are giants and subgiants with log^ < 3.5 as 
estimated in the KIC, so a reliable method for assigning ef- 
fective temperatures to such stars is highly desirable. Fortu- 
nately, this is feasible because the color-temperature relations 
for the bulk of the long cadence targets are not strong func- 
tions of surface gravity. For the purposes of the catalog we 
therefore supplement the fundamental dwarf scale with theo- 
retical corrections for the effect of surface gravity on the col- 
ors. 

Theoretical model atmospheres can be used to quantify the 
\ogg dependence of the color- temperature relations by com- 
paring the spectral energy distributions of dwarfs and giants. 
Figure |7] shows color-temperature relations along a 1 Gyr so- 
lar abundance isochrone for ^ - r, ^ - /, g- z, and J - Ks for 
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Fig. 7. — Main-sequence (solid; logg > 3.5) and post-main-sequence 
(daslied; logg < 3.5) color-temperature relationships for models along a 
1 Gyr isochrone with solar composition are shown for illustrative purposes. 
Colors illustrated are g - z, g - i, g — r, and J - Ks (from top to bottom). 



illustrative purposes. The model isochrone was taken from the 
web i nterface of the Padova isochrone database (Girardi et al. 
120021: iMarigo et al.l i2008fl As seen in Figure H] the model 
color- Teff relations are moderately dependent on logg, and il- 
lustrate that our photometric Jeff needs to be adjusted for gi- 
ants. 

We corrected for the difference in logg by taking theo- 
retical logg sensitivities in griz colors from the ATLAS9 
model atmosphere (Cast elli & Kuru cj 120041). The choice 
of these models seems internally inconsistent with our ba- 
sis model isochrone with MARCS-based colors. Neverthe- 
less, we adopted the ATLAS9 log g-color relations, primarily 
because our cluster-based empirical calibration of the color- 
Tgff relations has not been performed for subgiant and giant 
branches due to significant uncertainties in the underlying 
stellar interior models at these evolved stages. Therefore, it 
is just a matter of choice to adopt the ATLAS9 color tables in- 
stead of that of MARCS. Since we generated MARCS color 
tables in lAn et a l. (2009a) with a specific set of model param- 
eters for dwarfs (logg > 3.5), we simply opted to take the AT- 
LAS9 colors, and estimate a relative sensitivity of theoretical 
logg-color relations. 

We convolved synthetic spectra with the SDSS griz filter 
response curve fl and integrated flux with weights given by 
photon counts (Girardi et al.l2002h . Magnitudes were then put 
onto the AB magnitude system using a flat 363 1 Jy spectrum 
dOke & Gunnlll983l) . We created a table with synthetic colors 
from logg = 0.0 to 5.0 dex with a 0.5 dex increment, and from 
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Fig . 8 . — Theoretical Tiff corrections for various A log g values with respect 
to the fiducial isochrones. Corrections from Alogg = 0.5 to Alogg = 3.0 
with a 0.5 dex increment are shown. A linear ramp was used to define 
smoothly varying AT^ff over 4800 < T^ff < 5800 K. The sense is that giants 
with lower log g than the base isochi'one tend to have lower T^ff . 



4000 K to 6000 K with a 250 K increment at [M/H]=-1.0, 
-0.5, +0.0, and +0.2. Because YREC Tgff values were esti- 
mated at the fiducial metallicity, [Fe/H]= -0.2, we interpo- 
lated the color table to obtain synthetic colors at this metallic- 
ity. Note that Castelli & Kurucz (2004) adopted the solar mix- 
ture of Grevesse & Sauval (1998), as in ou r YREC isochrone 
model s (A09), so we assumed [M/H] in ICastelli & Kuruc j 
12004!) is the same as the [Fe/H] value. 

Figure [8] shows the correction factors in Teff as computed 
from synthetic colors as a function of colors in g - r, g - 
and g - z- We used our base isochrone to compute AT^ff at 
Alogg = 0.5, 1.0, 1.5, 2.0, 2.5, and 3.0 dex, where Alogg 
represents the difference between YREC logg and the logg 
in the KIC. The sense of ATeff is that giants with lower logg 
than the base model generally tend to have lower Teff than 
main-sequence dwarfs in the color range considered in this 
work. 

In Figure[8]we used a linear ramp over 4800 < Tgff < 5800 K 
(0.42 < g - r < 0.82), so that the theoretical ATeff becomes 
zero at T^ff > 5800 K. Otherwise the amplitude of theoretical 
Teff variations on the blue side {g - r < 0.6) would be similar 
to that of the red colors. Although this is not strictly true if the 
Alogg is large for blue stars, those stars are rare because stars 
on the giant branch (with the largest Alogg) have g - r> 0.5 
at near solar metallicity. The correction factors are tabulated 
in Table m If one wishes to adopt a different logg scale than 
in our base isochrone, tabulated ATeff factors can be used to 
correct for the logg difference. More importantly. Table |4] 
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TABLE 4 
\ogg Corrections 



Ref. Alogg 



g-r 


g-i/g-z 


logr' 


0.5 


1.0 


1.5 


2.0 


2.5 


3.0 


ATeff from g-r 


0.50 




4.57 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.55 




4.59 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.60 




4.61 


1.5 


2.3 


2.4 


2.0 


1.2 


-0.2 


0.65 




4.62 


5.8 


8.1 


7.3 


4.5 


-0.2 


-7.3 


0.70 




4.63 


11.2 


15.3 


12.9 


6.1 


-^.4 


-18.5 


0.75 




4.64 


18.5 


26.7 


24.6 


14.0 


-3.8 


-25.8 


0.80 




4.65 


28.2 


43.2 


44.3 


31.6 


7.3 


-22.8 


0.85 




4.65 


36.2 


58.6 


65.1 


55.1 


28.9 


-6.7 


0.90 




4.66 


41.2 


69.3 


81.1 


76.4 


54.4 


19.1 


0.95 




4.66 


44.6 


77.2 


93.7 


94.8 


78.5 


45.3 


1.00 




4.67 


47.1 


83.4 


103.8 


109.3 


97.9 


69.8 


1.05 




4.67 


50.6 


88.9 


110.8 


118.8 


111.4 


88.2 


1.10 




4.68 


51.5 


90.4 


111.8 


119.4 


114.0 


95.5 


1.15 




4.69 


53.1 


90.2 


106.4 


109.7 


104.0 


88.4 


1.20 




4.69 


48.4 


80.1 


88.0 


84.0 


76.0 


62.0 


AT^eff from g-i 


0.50 


0.655 


4.57 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.55 


0.725 


4.59 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.60 


0.795 


4.61 


1.2 


1.9 


1.8 


1.4 


0.6 


-0.5 


0.65 


0.865 


4.62 


5.2 


7.5 


6.6 


4.0 


-0.2 


-6.0 


0.70 


0.934 


4.63 


10.3 


14.5 


13.2 


8.2 


-0.4 


-12.0 


0.75 


1.003 


4.64 


16.9 


25.2 


25.8 


19.4 


7.0 


-10.3 


0.80 


1.071 


4.65 


24.7 


38.7 


43.4 


37.3 


22.4 


0.8 


0.85 


1.138 


4.65 


31.1 


50.8 


60.1 


57.1 


43.4 


21.5 


0.90 


1.205 


4.66 


34.9 


58.6 


71.4 


72.3 


61.9 


42.8 


0.95 


1.272 


4.66 


37.9 


65.1 


80.7 


85.2 


78.6 


63.0 


1.00 


1.341 


4.67 


41.3 


70.9 


88.3 


95.4 


92.1 


79.3 


1.05 


1.411 


4.67 


43.4 


74.3 


92.5 


101.1 


100.6 


90.4 


1.10 


1.483 


4.68 


43.8 


75.3 


93.6 


102.7 


103.9 


95.9 


1.15 


1.559 


4.69 


43.1 


73.7 


90.8 


98.7 


100.4 


94.2 


1.20 


1.639 


4.69 


38.3 


66.2 


80.6 


85.9 


86.7 


82.0 


Aleff fromg-z 


0.50 


0.708 


4.57 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.55 


0.795 


4.59 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.60 


0.881 


4.61 


0.8 


1.2 


1.1 


0.9 


0.4 


-0.4 


0.65 


0.966 


4.62 


3.6 


5.1 


4.4 


2.5 


-0.3 


-AA 


0.70 


1.050 


4.63 


7.7 


11.2 


10.1 


6.2 


0.5 


-7.5 


0.75 


1.133 


4.64 


13.5 


20.6 


20.7 


15.6 


7.5 


-4.2 


0.80 


1.215 


4.65 


20.4 


32.5 


35.5 


30.6 


21.0 


6.4 


0.85 


1.294 


4.65 


25.5 


41.7 


48.3 


46.5 


37.4 


22.5 


0.90 


1.373 


4.66 


28.1 


46.9 


56.5 


58.3 


50.9 


37.4 


0.95 


1.452 


4.66 


30.2 


51.5 


63.8 


68.4 


63.3 


51.8 


1.00 


1.533 


4.67 


33.0 


56.5 


70.5 


76.7 


74.0 


64.3 


1.05 


1.616 


4.67 


34.8 


59.9 


74.9 


82.2 


81.7 


73.9 


1.10 


1.703 


4.68 


35.4 


61.3 


76.7 


84.6 


85.9 


80.0 


1.15 


1.794 


4.69 


34.6 


60.0 


75.2 


83.1 


85.4 


81.2 


1.20 


1.891 


4.69 


30.3 


53.7 


67.8 


74.9 


77.6 


75.2 



Note. — The sense of the difference is that a positive ATefi- means a higher r^fr at a lower logg. 
^ The logg values in the YREC model. 



can be used to infer Tgff for giants, since our base isochrone 
(Table [Til covers stellar parameters for main-sequence dwarfs 
only. 

The biggest AT^s in Figure [8] is ^ 100 K. However, the 
effects of the log^ corrections are moderate in the KIC. If 
we take the mean ATeff correction in g - r, g - i, and g - z, 
the mean difference in T^ff between KIC and YREC decreases 
from 190 K to 166 K for stars with logg < 3.5. The logg 
corrections are insensitive to metallicity. The ATeff in Fig- 
ure[8]was computed at [Fe/H]= -0.2, but these corrections are 
within 10 K away from those computed at [Fe/H]= -0.5 {^la 
lower bound for the KIC sample) when A logg = 1. 

The statistical properties of the SDSS giant temperatures 
are compared with spectroscopic data in Section [J!4l and with 
the KIC in Section ffl 



3.3. Tests with Open Cluster Data 

The IRFM technique provides global color-metallicity-rgff 
correlations using field samples, while clusters give snapshots 
at fixed composition, which define color- Teft trends more pre- 
cisely. Deviations from color to color yield the internal sys- 
tematic within t he system , as the color-temperature relation- 
ships defined in lAn et al.l (l2007bl) are empirical descriptions 
of actual cluster data. The A09 SD SS system, by construction, 
agrees with the lAn et al.l (12007 b') Johnson-Cousins system; 
but we can check the concordance between the two scales 
within the open cluster system. 

We have two basic results from this comparison. First, 
(J - Ks)-hased temperatures from the IRFM are different from 
other IRFM thermometers. J - Ks is also the only IRFM di- 
agnostic available for the bulk of the KIC sample. When ac- 
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Fig. 9. — Internal consistency of the IRFM T^ff estimates for the Hyades 
(red circles) and Praesepe stars (blue triangles). Comparisons are shown for 
each color index with respect to the T^ff values determined from V — /c at 
[Fe/H]= 0.13 for the Hyades and [Fe/H]= 0.14 for Praesepe. En'or bars rep- 
resent itlcr uncertainty propagated from photometiic en'ors. 

counting for the offset in J -Kg relative to other IRFM indica- 
tors, the underlying IRFM system and the SDSS system are in 
excellent agreement for stars below 6000 K. Second, there is 
a systematic offset between the IRFM and SDSS scales above 
6000 K. We therefore correct the high end temperature esti- 
mates for the SDSS to put them on the IRFM scale, which 
yields an internally consistent set of photometric temperature 
estimates. 

Figures |9HTT] show how the IRFM Teff determinations are 
internally consistent in the Jounson-Cousins-2MASS system 
in B - V, V — Ic, V - Ks, and J - using stars in four well- 
studied clusters: The Hyades (red circles in Figure |9l), Prae- 
sepe (blue triangles in Figure |9l), the Pleiades (Figure [TOl l. 
and M67 (Figure [TTl i. All of the stars shown in these fig- 
ures are likely single-star members of each cluster after ex- 
cluding known (unresolved) binaries. In FigurefTT] we show 
results based on the two independent sets of M67 photometry 
from Montgomery et al. ( 1993, blue triangles) and SandquisJ 
(l2()04, red circles). The compilation an d individu al sources of 
the cluster photometry can be found in I An et aP f2007b ). 

To construct Figures l9HTn we corrected observed mag- 
nitudes for extinction using E{V - Ic)/E(B - V) = 1.26, 
E(V - K,)/E(B - V) = 2.82, and E(J - K,)/E(B - V) = 
0.53 ( An et al.l l2007al) . Foreground reddening values of 
E(B -V) = 0.000 ± 0.002, 0.006 ± 0.002, 0.032 ± 0.003, 
and 0.041 ±0.004 mag were used fo r the Hyades, Pra esepe, 
the Pleiades, and M67 respectively (l An et al.ll2007bl) . The 
IRFM Teff equations in CIO include metallicity terms, and we 
adopted [Fe/H]= -1-0.13 ± 0.01, -l-0.14±0.02, -1-0.04 ± 0.02, 
and -1-0.00 ± 0.01 dex for the Hyades, Praesepe, the Pleiades, 
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Fig. 10. — Same as in Figure m but for the Pleiades at [Fe/H]= 0.04. Note 
that low-mass Pleiades stars (T^ff < 5000 K) are known to have anomalously 
blue colors in B — V. These stars could also have slight near-IR excesses, 
which may have affected Teff values from J — Ks. 
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and M67, respectively, based on high-re solution spectro - 
scopic abundance analysis (see references in lAn et al.l2007bl) . 
Only the (B - y)-based estimates are significantly impacted 
by metallicity corrections, and the relative abundance differ- 
ences in these well-studied open clusters are unlikely to be 
substantial enough to affect our results. 

The ±1ct error bars in Figures IgHTTI are those propagated 
from the photometric errors only. Mean differences in the 
IRFM Teff and the errors in the mean are provided in Ta- 
ble|5] Global differences are shown for stars at 4000 < Teff < 
7400 K, and those cooler and hotter than 6000 K are shown 
in the table. The CTsys represents a total systematic error in 
this comparison from the reddening and metallicity errors 
(summed in quadrature); however, systematic errors are less 
important than random errors because of the precise E(B - V) 
and [Fe/H] estimates of these well-studied clusters. 

The low-mass stars in the Pleiades are known to have 
anomalously blue colors related to stellar act ivity in these 
heavi ly spotted, rapidly rotating, young stars (IStauffer et al.l 
I2003h . The temperature anomaly for B - V at Teff < 5000 K 
in Figure [TOl which is ^ 200 K larger than that for more 
massive stars, reflects this known effect and therefore is not 
a proper test of internal consistency in old field stars (such 
as those in the KIC). The M67 data may also be inappropri- 
ate for the test of the IRFM internal consistency, but with a 
different reason. Two independent photometry sets lead to 
a different conclusion: iMontgomerv et al] {1993) photometry 
shows inte rnally less consistent IRFM Teff for M67 stars than 
[Sandauist (2004). A similar argument was made in An et al.l 
(l2007b), based on the differential metallicity sensitivities of 
stellar isochrones in differe nt color indices (se e Figure 11 in 
the above paper); see also IVandenBerg et al.l (l2010i) for an 
indepen dent confirmation of the s ystematic zero-point issue 
with the Montg omery et al.l (Il993h photometry. 

Our cluster tests based on the Hyades and Praesepe demon- 
strate the internal consistency of the CIO color- Teff relations 
in B - V, V - Ic, and V - K,. The mean differences in T^ff 
among these color indices are typically few tens of degrees 
for both hot and cool stars (Table |5]l. Some of these mean 
differences could be systematic in nature, but they are gener- 
ally consistent with the scatter in the CIO IRFM calibrations. 
However, the (J - /T,)- Teff relation tends to produce hotter 
Teff than those from other color indices for these cluster stars 
(see bottom panel in Figure|9]l. The mean differences between 
Teff(V - Ic) and 7;ff(7 - K,) are 95 K and 44 K for the Hyades 
and Praesepe, respectively. There is also a hint of the down- 
turn in the comparison for the hot stars in these clusters, where 
(J - Ks)- Teff produces cooler temperatures than {V - Ic) - Teff 
relation. The 100 K offset between the hot and the cool 
stars roughly defines the size of the systematic error in the 
IRFM technique of CIO in 7 - 

The Pleiades stars show a weaker systematic Teff trend for 
the cool and the hot stars than the Hyades and Praesepe. In 
spite of this good agreement, we caution that this could be a 
lucky coincidence because the Pleiades low-mass stars prob - 
ably have slight near-lR excesses in K, (IStauffer et al.ll2003h . 
The main-sequence turn-off of M67 is relatively cool, so the 
difference is only suggestive. 

Figure [12] shows comparisons between the IRFM and 
YREC Teff estimates. The left panel shows the comparisons 
for the Hyades and Praesepe stars, while the right panel shows 
those for the Pleiades and M67 stars. The IRFM Teff on top 
panels was computed based on the (V - Ic) - 7"eif relation in 
CIO, just as those used for a principal Teff estimator in the 



above compar isons (Fig ures IgfflHi. The YREC T^^ff was es- 
timated using lAn et al.l (l2007b) isochrones, which have the 
same underlying set of interior models as those used in the 
current analysis. The model T^ff was computed at a con- 
stantMy of individual stars, assuming (m - M)o = 3.33 ±0.01, 
6.33 ±0.04, 5.63 ±0.02, and 9.61 ±0.03 mag for the dis- 
tance moduli of the Hyades (550 Myr), Praesepe (550 Myr), 
the Pleiades (100 Myr), and M 67 (3.5 Gyr), respectively (see 
references in An et al.ll2007bh . 

Table |5] lists weighted mean differences between YREC 
and IRFM Teff. The mean difference between the {V - Ic)- 
based IRFM and the luminosity-based YREC Teff for cool 
stars (Teff < 6000 K) is less than 20 K, but the differences 
rise above 6000 K to the 50 K level. The difference between 
the (J - /:,)-based IRFM and M^-based YREC T^ff shows dif- 
ferent offsets for the cool and hot stars; this trend is consistent 
with the above comparison between (J - Ks)-hased IRFM and 
other IRFM determination. 

The bottom panel in Figure[T2lshows comparisons between 
the YREC T^n and the average IRFM T^n from B-V,V -Ic, 
and V - K,. Our results using J - as a thermometer are 
consistent with our earlier finding in Section [TTI that CIO 
(J - Ks)-hased T^ff values are systematically cooler than those 
from the griz-hased YREC models for hot stars (above about 
6000 K). The (J - K,)-hased differ both from other IRFM 
diagnostics and the values inferred from SDSS colors for 
cooler stars, while the mean values inferred from the IRFM 
are close to SDSS for the cooler stars. 

We therefore conclude that the cool star temperature scales 
are consistent, while there is evidence for a systematic de- 
parture at the hot end. A similar pattern emerges when we 
compare with spectroscopy, as discussed in the next section. 
Caution is therefore required in assigning errors for stars with 
formal temperature estimates above 6000 K. 

Systematic Tgff differences are shown in Figure [13] The red 
line represents the difference with the { J - ^r5)-based IRFM 
Teft for the open cluster sample (Hyades and Praesepe), while 
the orange line shows that with respect to the mean IRFM val- 
ues from B - V ,V - Ic,V - K^. Error bars indicate ±1(T error 
in the mean difference. The difference between the average 
IRFM scale and the SDSS scale in the clusters is less than 
25 K on average from 4000-6000 K, which we take as a con- 
servative systematic temperature uncertainty in that domain. 
The differences are moderately larger for the IRFM J - K, 
temperature alone, but that diagnostic is also different from 
other IRFM thermometers for cool stars. 

The differences in the hot cluster stars reflect actual differ- 
ences in the calibrations, not issues peculiar to the photome- 
try, extinction, or blending. We therefore attribute the compa- 
rable differences seen in the KIC stars (gray band) as caused 
by calibration issues in J - Ks rather than as a reflection of 
systematics between the IRFM and SDSS systems. Further- 
more, the SDSS calibration was based on M67 data, where 
the hotter turnoff stars (T^ff > 6000 K) were saturated. As a 
result, we believe that an adjustment closer to the IRFM scale 
is better justified. 

A simple correction term, of the form below 

Teff^sDSS < 6000 K : nif,,,„ = 7;ff(SDSS), (6) 
6000 K < reft;sDSS < 7000 K : (7) 
Teff.corr = 0.8 (reff(SDSS) - 6000 K) + 6000 K, (8) 
Teff.SDSS > 7000 K : Teff.con = TefKSDSS) - 200 K (9) 
brings the two scales into close agreement across their mu- 
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TABLE 5 

Statistical Properties of Clusters Comparisons 



Afeff (K) 



Cluster Data 


4000 < T^fi < 7400 




4000 < T^ff < 6000 


6000 < Vrff < 7400 






Teff(B - 


V,IRFM)- 


Teff(V - 


/c.IRFM) 






Hyades 


10.3 ±6.7 






11.5±6.0 


-2.9 ±15.3 


4.0 


Praesepe 


11.3±7.4 






23.0 ±7.4 


-59.4±15.3 


6.7 


Pleiades 


54.3 ±12.1 






67.5 ±13.3 


-5.5 ±19.0 


6.5 


M67 (MMJ93)'' 


74.1 ±7.9 






71.5±8.5 


93.2±21.0 


4.3 


M67 (S04)'' 


28.8±4.0 






16.7 ±4.2 


36.1 ±7.7 


5.3 




Trff(y - 


^Cs.IRFM)- 


Tcff(V - 


/c,IRFM) 






Hyades 


-3.8 ±7.4 






-4.0 ±7.6 


-1.9±16.7 


0.9 


Praesepe 


23.2 ±6.2 






25.0 ±6.7 


1.6±13.7 


1.0 


Pleiades 


10.8±5.9 






9.0 ±6.9 


19.2±13.9 


1.1 


M67 (MMJ93)'' 


87.3 ±7.7 






81.8±8.7 


133. 8± 19.6 


0.5 


M67 (S04)'' 


34.1 ±3.9 






36.7 ±4.9 


29.9 ±6.2 


1.2 




Teff(y - 


a:,,irfm)- 


Teff(V - 


/c.IRFM) 






Hyades 


94.7 ± 17.0 






108.3 ±15.7 


36.0 ±32.9 


3.9 


Praesepe 


44.2 ±13.1 






66.2 ±13.7 


-97.8 ±34.8 


4.0 


Pleiades 


-10.1 ±13.9 






-0.6 ± 14.8 


-55.3 ±30.7 


5.4 


M67 (MMJ93)'' 


51.8±16.6 






61.4±18.3 


9.1 ±38.6 


6.3 


M67 (S04)'' 


-47.3 ±14.1 






-18.7± 18.5 


-86.7±21.7 


7.3 




Teff(V 


-/c.IRFM) 


-Teff(M 


v,YREC) 






Hyades 


-12.1 ±8.7 






-8.3 ±8.3 


-56.4 ±15.4 


6.4 


Praesepe 


-14.7±6.6 






-12.2 ±7.4 


-52.9 ±11.9 


21.3 


Pleiades 


-14.6 ±8.4 






-9.9 ±8.4 


-50.5 ±20.3 


13.9 


M67 (MMJ93)'' 


-52.9 ±8.6 






-34.4 ±7.5 


-178.4±13.7 


22.0 


M67 (S04)'' 


-20.4 ±3.4 






-13.5 ±3.5 


-37.1 ±8.2 


24.1 




Teff(7- 


-/f,,IRFM) 


-Teff(Mv/,YREC) 






Hyades 


62.2 ± 16.6 






90.9 ±12.5 


-35.1 ±32.5 


6.9 


Praesepe 


7.9 ± 14.9 






43.9 ±13.0 


-152.5 ±27.5 


17.5 


Pleiades 


1.2±15.9 






22.6±17.6 


-92.8 ±26.9 


11.4 


M67 (MMJ93)'' 


-50.1 ±15.9 






-16.5 ±19.2 


-123.8±28.4 


35.3 


M67 (S04)'' 


-60.1 ±15.8 






-31.9 ±18.8 


-127.7±29.1 


30.4 



^ Systematic errors from reddening and metallicity, summed in quadrature. In the comparisons between IRFM and YREC, we also include effects of the cluster age and distance 
modulus enors. 
MMJ= :Montgomerv et aljfT993l) ^ S04 jSandquisiHSoM) . 
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Fig. 12. — Left: The Teff comparisons between IRFM and YREC for stars in the Hyades (red circles) and Praesepe (blue triangles). Right: Same as in the 
left panels, but for the Pleiades stars (black circles), M67 stars from Sandquist 1 2004, red circles) photometry and from that of Montgomery et al. 1 1993, blue 
triangles). Top: IRFM T^ff from V - Ic. Bottom: Mean IRFM Teff from B-V,V -Ic, and V - K,. The YREC T^fi was estimated from the luminosity (Mv) of 
each star. Black line in the bottom left panel shows a moving averaged trend of the T^ff difference. 
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Fig. 13. — Sy.stematic difference.s of various T^n estimates with respect to tlie YREC scale. Grey line shows the mean trend for the main KIC sample discussed 
in this work. The red line represents the difference with the (J — ^r5)-based IRFM Teff for the open cluster sample (Hyades and Praesepe), while the orange line 
shows that with respect to the mean IRFM values from B -V.V — Ic,V — Kg. The blue line shows the trend for the Bl 1 spectroscopic sample. Error bars in 
all cases represent itlo" en'or in the mean difference. Our adopted hot-Teff coiTections are shown with a black dashed line. Note that the empirical color- Teff 
conections in YREC are defined at 4000 < T^ff < 6000 K in SDSS colors. 
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tual range of validity. This empirical correction is indicated 
by the black dashed line in Figure [T3] Below we find offsets 
similar in magnitude and opposite in sign between the IRFM 
and spectroscopic temperatures for hotter stars. Although this 
does not necessarily indicate problems with the fundamental 
scales, it does imply that systematic temperature scale differ- 
ences are important for these stars. 

3.4. Comparison with Spectroscopy 

Spectroscopy provides a powerful external check on the 
precision of photometric temperature estimates. Spectro- 
scopic temperatures are independent of extinction, and can 
be less sensitive to unresolved binary companions and crowd- 
ing. In this section we therefore compare the photometric 
and spectroscopic temperature estimates for t wo w ell-studied 
samples in the Kepler fields. iBruntt et al. I (1201 1 , hereafter 
Bll) reported results for 93 stars with asteroseismic data, 
includ ing 83 stars in our sample. Mo lenda-Zakowicz et alj 
(1201 IL hereafter MZll) reported results for 78 stars, includ- 
ing 45 targets in common with our sample. The MZll data 
for cool stars are mostly subgiants and giants, while the bulk 
of the dwarf sample is hotter than 6000 K. The Bll sample is 
similarly distributed, with the transition from the cool evolved 
to the hot unevolved sample occurring at 5500 K. 

All comparisons below are for the corrected photometric 
scale, adjusted for concordance with the IRFM at the hot 
end. We compare spectroscopic methods both with the fixed- 
metallicity ([Fe/H]= -0.2) temperatures in the catalog and the 
refined temperature estimates made possible with the addition 
of metallicity information and theoretical metallicity correc- 
tions. We excluded outliers in the following statistical com- 
parisons using a 3cr outlier rejection. 

As demonstrated below, we find that the two spectroscopic 
samples have different zero-points with respect to both the 
SDSS and KIC samples, indicating the importance of system- 
atic errors in such comparisons. The photometric scale for the 
cool dwarfs and giants are in good agreement with the Bll 
scale, while both are offset relative to MZll. The situation 
is different for hot dwarfs. The IRFM scale was cooler than 
the uncorrected SDSS scale. The spectroscopic samples are 
cooler than both. We interpret this as evidence of additional 
systematic uncertainties for the F stars, and discuss possible 
causes. 

The stellar parameters for the MZl l samp le were de- 
rived using the ' Molenda-Zakowicz et akl (1200 7*) template ap- 
proach. The spectra were compared with a library of reference 
stars. The surface gravity, effective temperature, and metallic- 
ity were derived from a weighted average of the five closest 
spectral matches in the catalog. Bll used asteroseismic sur- 
face gravities and derived effective temperatures from tradi- 
tional Boltzmann-Saha consistency arguments. 

We compare the spectroscopic and photometric tempera- 
ture estimates in Figure[T4l The top, middle, and bottom pan- 
els compare spectroscopic temperatures to those of the KIC, 
IRFM (7 - Ks), and SDSS, respectively. Left panels show 
comparisons for dwarfs (KIC logg > 3.5), while the right pan- 
els show those for giants (KIC \ogg < 3.5). Filled circles are 
the Bll data, while open circles are the MZl 1 data. In total, 
83 out of 93 sample stars in Bll were used in this compari- 
son; the remaining 10 stars do not have griz photometry in all 
passbands, so were not included in our KIC subsample. For 
the same reason, we initially included 45 spectroscopic targets 
from MZl 1 , but later excluded 8 more stars with g - r < 0. 1 or 



g — r> 1.0. Triangles in the bottom two panels represent stars 
flagged as having internally inconsistent effective temperature 
estimates (Section [4.2.41 ). Error bars show the expected ran- 
dom errors, with a 70 K error adopted in the temperature for 
the individual Bll sample stars. 

In the above comparisons, we corrected the IRFM temper- 
ature estimates for the spectroscopic metallicity measurement 
of each sample, although the T^ff corrections in CIO were 
negligible (ATeff « 18 K) in J - Kg. We also used individ- 
ual stellar isochrones at each spectroscopic metallicity to es- 
timate SDSS Teff from griz, assuming a constant age of 1 Gyr 
at all metallicity bins. However, the net effect of these cor- 
rections was small (ATss w 25 K), because griz-T^s relations 
are relatively insensitive to metallicity and the mean metal- 
licities of the spectroscopic samples are close to our fiducial 
value (([Fe/H]) = -0.07 ± 0.02 and -0.1 1 ± 0.03 for the Bl 1 
and MZll samples, respectively). The SDSS T^ff values for 
giants were corrected for the \ogg difference from the dwarf 
temperature scale as described in Section [372l 

Both spectroscopic samples for dwarfs are systematically 
hotter than the KIC (top left panel in Figure[T4li. The weighted 
average difference between the Bll sample and the KIC, in 
the sense of the KIC minus spectroscopic values, is -170 K 
with a dispersion of 116 K, after a 3o- outlier rejection. The 
MZl 1 sample is closer to the KIC, with a -82 K mean dif- 
ference and a dispersion of 172 K. This difference of 88 K is 
a reflection of the systematic errors in the spectroscopic tem- 
perature scales. In the above comparisons, we did not include 
stars with inconsistent SDSS temperature measurements (tri- 
angles in Figure O. 

The weighted average difference between the Bll sample 
and the SDSS (in the sense SDSS - Spec) for dwarfs is 85 K 
with a 95 K dispersion, after excluding those flagged as hav- 
ing discrepant reff(YREC) values. If the metallicity correc- 
tions to the SDSS values were not taken into account (i.e., 
based on models at [Fe/H]= -0.2), the mean difference be- 
comes 73 K, but the dispersion increases to 1 1 1 K. 

However, there is a strong temperature dependence in the 
offset. Below 6000 K the mean difference is 50 K with a 
dispersion of 47 K. For the hotter stars the mean difference 
is 101 K with a dispersion of 118 K. The blue line in Fig- 
ure [T3] shows a moving averaged difference between the Bll 
spectroscopic values and SDSS T^n without the hot-end Teff 
corrections (equations 6-8). 

Although the size of the dwarf sample in MZll is small, 
it is found that the effective temperatures are systematically 
cooler than the SDSS values, with a weighted mean offset of 
152 K (SDSS - Spec) and a dispersion of 175 K. The dif- 
ference is temperature dependent, being 53 K for the stars 
below 6000 K and 178 K above it. These differences are 
3 K and 77 K larger, respectively, than the results from the 
Bll sample. The temperature differences between photome- 
try and spectroscopy are therefore smaller than the differences 
between the spectroscopic measurements and the KIC, while 
there is a real difference at the hot end even when system- 
atic differences between the two spectroscopic samples are 
accounted for. 

The Bl 1 sample includes only two giants (KIC logg < 3.5), 
but their spectroscopic temperatures are consistent with both 
IRFM and SDSS temperatures (see middle and bottom right 
panels in Figure O. On the other hand, the MZll sample 
shows a large offset from IRFM (A7;i=f = 245 K) and SDSS 
( A^eff = 206 K), while the KIC and the MZl 1 values agree 
with each other (A7;ff = 9 K). 
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giant comparisons (KIC logg < 3.5). Triangles in the bottom two panels represent stars flagged as having internally inconsistent effective temperature estimates 
(see text). 



The cool MZl 1 stars are mostly subgiants and giants, while 
the B 1 1 cool sample includes a large dwarf population be- 
tween 5000 K and 6000 K. The difference between the two 
cool end results - good agreement with Bll for cool dwarfs, 
but not with MZll - is real. This could reflect systematic 
differences between the dwarf and corrected giant results for 
the SDSS or the templates adopted by MZl 1 for the evolved 
and unevolved stars. The scatter between the MZll results 
and the photometric ones is substantially larger than that be- 
tween Bll and photometric temperature estimates. It would 
be worth investigating the zero-point of the templates used in 
the former method, as well as the random errors, in light of 
the results reported here. 

In the section above we have focused on differences be- 
tween the scales; it is fair to ask how both might compare to 
the true temperatures. The photometric scale is at heart sim- 
ply an empmcal relationship between color and the definition 
of the effective temperature itself (L = AnR^aT^f^, and there- 
fore the scale itself should be sound where the photometric 
relations are well-defined. However, the photometric methods 



can fail if there is more than one contributor to the photome- 
try, or if the reddening is incorrectly measured. Spectroscopic 
temperatures measure physical conditions in the atmosphere, 
and are only indirectly tied to the fundamental flux per unit 
area, which defines the effective temperature. There are also 
systematic uncertainties between different methods for infer- 
ring effective temperatures, for example, fitting the wings of 
strong lines, or the use of Boltzmann-Saha solutions based on 
ionization and excitation balance. Finally, both photometric 
and spectroscopic estimates are only as good as their assump- 
tions; stars with large surface temperature differences will be 
poorly modeled by both methods. 

Our primary conclusion is therefore that the various dwarf 
temperature methods, spectroscopic and photometric, are in 
good agreement for the cooler stars. Systematic effects are 
at or below the 50 K level. The hotter stars in the sample 
have real systematic differences between spectroscopic and 
photometric temperatures, and similar discrepancies are also 
present between the photometric methods themselves. This is 
further evidence that work is needed to tie down more pre- 
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cisely the temperature scale above 6000 K, and that larger 
systematic errors should be assigned in this domain until such 
an analysis is performed. We have less data for the giants, 
but there does appear to be a real difference between the pho- 
tometric results and the temperatures inferred for the MZl 1 
sample. 

3.5. Effects of Binaries on Colors 

Unresolved binaries in the sample could bias a color-based 
Teff estimate. Unless the mass ratio of the primary and 
secondary components in the binary system is close to ei- 
ther unity (twins) or zero (negligible contributions from the 
secondary), composite colors of the system are redder than 
those from the primaries alone, leading towards systemati- 
cally lower Tgff. It is difficult to directly flag potential bina- 
ries given the filters available to us, and as a result we do not 
include star by star corrections in the table. However, such 
a systematic bias will be important when evaluating the bulk 
properties of the KIC sample. In this section, we therefore 
estimate the size of the bias due to unresolved binaries in the 
KIC, and provide statistical corrections for the effect of un- 
resolved binary companions on average effective temperature 
estimates. 

Binary contamination effects on the color- Teff relations 
were derived by performing artificial star tests . We used a 
1 Gy r old Padova models at solar abundance dGirardi et alj 
120041) . These models include stellar masses down to 0.15 Mq, 
allowing us to include low-mass systems outside the formal 
range of the SDSS color calibration. The absolute color- 
Teff relations in these models are not exactly the same as in 
our base calibration, and the adoption of a solar metallicity 
isochrone is not strictly self-consistent with our application 
of the base model at [Fe/H]= -0.2. However, our main pur- 
pose is to evaluate the relative temperature errors induced by 
companions, and the effects of these offsets are presumably 
small. 

We assumed a 50% binary fraction with 10,000 single stars 
and 10,000 binary systems. Primary masses were randomly 
drawn from a Salpeter mass function, while we explored 
three different choices for the relative masses of the secon- 
darie s: Salpeter, flat, and one drawn from the open cluster 
M35 (iBarrado v Navascues et al.ll2001h . A flat mass function 
is expected for short-period binaries, which will be a minor- 
ity of the sample; this is thus a limiting case. In the artificial 
star simulations, we derived empirical color-color sequences 
in g - i, g - z, and J - with g - r as the principal color 
index. We simulated photometric errors by injecting disper- 
sions of 0.01 mag in gri, 0.03 mag in z, 0.024 mag in J, and 
0.028 mag in K^. These 2MASS errors are median values of 
the actual photometric errors in the KIC sample. 

The result of these binary simulations is presented in Fig- 
ure [15] which shows the mean deviations in g - i, g - z, and 
J - Ks from those with primaries alone. For Figure[T5]we fit- 
ted a Gaussian for each g - r bin to estimate the mean color 
offset and the uncertainty as shown by circles and error bars. 
The three curves indicate results from three different relative 
mass functions for secondaries. 

The sizes of these color shifts are shown in Table |6l 
The systematic color shift due to unresolved binaries is less 
strongly dependent on the choice of secondary mass func- 
tions. Typical sizes of these color shifts are ^ 0.003 mag, 
0.008 mag, and 0.010 mag in ^ - /, and g- z, and J - K,, re- 
spectively. To correct for the unresolved binaries in the KIC, 
the above color shifts should be subtracted before estimating 
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Fig. 1 5 . — Average color bias in g — i, g — z, and J - Ks at fixed g — r due to 
unresolved binaries for three different assumptions about the secondary mass 
function. Points and error bars are the centroid and the error in the mean 
distribution from the simulations. A 50% binary fraction is assumed. 

Teff. The last four columns in Table |6] list the average Teff dif- 
ference between a population with a 50% unresolved binary 
fraction and that of primaries alone. The sense is that un- 
resolved binary stars have lower temperatures than expected 
from primaries alone. Different SDSS color indices have sim- 
ilar binary sensitivities, and temperatures based on these fil- 
ters are less affected by unrecognized companions than those 
derived using J - Ks- These color shifts are small for any 
given star, but significant when applied to the entire catalog. 
We therefore recommend including them when using large 
samples of photometric effective temperature estimates, and 
include this effect in our global error budget below. 

3.6. Other Sources of Uncertainties and Error Budget 

We can assess our overall errors by comparing the real to 
the observed dispersions in the color-color plane. Photomet- 
ric errors, unresolved binaries, and metallicity all induce scat- 
ter; so would extinction uncertainties. Significant mismatches 
between the two reflect unrecognized or overestimated error 
sources. 

Figure [16] shows the observed color-color diagrams in the 
KIC, after the extinction corrections and the zero-point ad- 
justment as described in Section 12.21 From Figure [16] we 
estimated the standard deviation of the color dispersion from 
a fiducial line (fit using a 5''' order polynomial) in g - i, g - z, 
and J - Ks at each g - r bin. These observed dispersions 
with good Teff estimates are shown as solid black curves with 
closed circles in Figure [17] Here the criteria for the good Teff 
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TABLE 6 
Binary Corrections 





A(g - i) 


A(g-z) 




g-r 




g-z 


J-K, 


(mag) 


(mag) 


(mag) 


(mag) 


(K) 


(K) 


(K) 


(K) 


M35 Mass Function'' 


0.25 


0.001 


0.000 


0.006 


18 


21 


23 


83 


0.35 


0.000 


0.002 


0.012 


25 


28 


32 


89 


0.45 


0.002 


0.004 


0.010 


27 


32 


35 


70 


0.55 


0.004 


0.006 


0.010 


31 


38 


43 


77 


0.65 


0.002 


0.007 


0.009 


32 


36 


42 


93 


0.75 


0.002 


0.010 


0.014 


30 


34 


40 


121 


0.85 


0.002 


0.012 


0.008 
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37 
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Note. — The sense of the bias is that populations mixed with unresolved binaries look redder (cooler) at a given ^ — r in the above color indices. 
^ Mean difference in Tefr between a population with a 50% unresolved binary fraction and that of primaries alone. The sense is that unresolved binary stars have lower temperatures 
than expected from primaries alone. 

" Mass function for secondary components in the binary system. All simulation results are based on a 50% unresolved binaiy fraction. 



are that the standard deviation of individual Tgff from three 
color indices {g- r, g- i, and g - z)is less than 130 K or that 
the difference between SDSS and IRFM measurements is no 
larger than three times the ra ndom errors of these measure- 
ments (see also Section |4. 2. 4l i. There is a strong overlap be- 
tween the two criteria. Since the formal random SDSS errors 
are of order 40 K, and the systematics between the colors are 
typically at that level as well, differences of 130 K represent 
clear evidence of a breakdown in the color-temperature rela- 
tionships, likely from unresolved blends. Excluding extreme 
outliers is essential because they would otherwise dominate 
the dispersion measure, and we are interested in testing the 
properties of the majority of the sample. 

Other lines in Figure [17] represent contributions from ran- 
dom photometric errors, unresolved binaries or photometric 
blends, metallicity, and dust extinction as described below. 
Red lines with open circles are the quadrature sum of all of 
these error sources. 

We assumed 0.01 mag errors in gri, 0.03 mag errors in 
z, 0.024 mag in /, and 0.028 mag in K, to estimate color 
dispersions from photometric errors alone (red curve in Fig- 
ure [iTll. To perform this simulation in grizJKs, we combined 
our base model (Table [D with our earlier set of isochrones 
in the 2MASS system (An et al. 2007b|3 at the same metal- 
licity ([Fe/H]= -0.2) and age (1 Gyr) as those for the base 

'" Available at |http://www.as tronomy.ohio-state.edu/iso/pl.html | 



isochrone. As with the binary simulations described in the 
previous sectio n, we employed a 1 Gyr old, solar metallicity 
Padova model dGirardi et al.ll2004l) to generate color-color se- 
quences with a 50% binary fraction based on the M35 mass 
function for secondaries. Again, running this isochrone in the 
simulation is not strictly consistent with the usage of our base 
model, but the relative effects induced by unresolved com- 
panions would be rather insensitive to the small metallicity 
difference. The dispersion induced by unresolved binaries is 
shown in a red dashed curve in Figure [TT] 

The KIC sample has a mean [Fe/H]= -0.2 with a standard 
deviation of 0.28 dex. If the KIC [Fe/H] values are accu- 
rate enough for these stars, this metallicity spread would in- 
duce a significant spread in T^ff. The color dispersion due to 
metal abundances was estimated by taking the color differ- 
ence between our base model ([Fe/H]= -0.2) and the mod- 
els at [Fe/H]= +0.1 and -0.5 as an effective ilcr uncer- 
tainty. The metallicity error contribution is shown in blue 
solid curves. The KIC sample has a wide range of redden- 
ing values (0 < E(B -V) < 0.2). We took 0.02 mag error 
as an approximate ilcr error in E{B - V), roughly equivalent 
to a 15 % fractional uncertainty for a typical star. Stars on 
the simulated color-color sequence were randomly displaced 
from their original positions assuming this E(B - V) disper- 
sion. The resulting color dispersion is shown with the blue 
dashed curves in Figure [TT] 

In Figure [17] there is a color-dependent trend in the error 
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Fig. 16. — Extinction corrected color-color relat ions in the KIC, after the 
zero-point corrections as described in Section [T2l Only those with logg > 
3.5 are shown. 



budget, where observed color dispersion increases for cooler 
stars in g - / and g- z- On the other hand, the simulated dis- 
persions (open circles connected with solid red curves) are 
essentially flat. Our results are consistent with expectations 
m J - K; if anything, the random errors appear to be over- 
stated. This is probably caused by correlated errors in J and 
Ks, which were treated as uncorrected in the temperature er- 
ror estimates. 
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Fig. 17. — Comparison between observed (thick black line with closed cir- 
cles) and modeled (thick red line with open circles) dispersions of the color- 
color sequence as a function of g — r. The modeled dispersion is a quadrature 
sum of individual error contributions: photometric errors (red solid), unre- 
solved binaries (red dashed), metallicity (blue solid), and reddening (blue 
dashed). 

Based on this exercise, we conclude that our error model is 
reasonable for the hot stars in the sample, especially when the 
stars most impacted by blends are removed. There is excess 
color scatter for red stars, which correspond to effective tem- 
peratures below ^ 5000 K in our sample. About 16% of the 
sample are found in this temperature domain. This could re- 
flect contamination of the dwarf sample by giants, which have 
different color-color relationships; or a breakdown in the pho- 
tometric error model for red stars. It would be useful to revisit 
this question when we have a solid estimate of the giant con- 
tamination fraction for the cool dwarfs in the sample. 

4. THE REVISED Teff CATALOG 

4.1. A Recipe for Estimating Teff 

We present results for the long-cadence sample with the 
overall properties of the catalog and systematic error esti- 
mates in this section. We have not provided corrected values 
for the entire KIC, because the additional quality control is 
outside the scope of our effort. However, our method could be 
applied in general to the KIC, employing the following steps. 

• Correct the KIC griz photometry onto the SDSS DR8 
system using equations 1-4. 

• Apply the KIC extinctions and the extinction coeffi- 
cients in Section|2]to obtain dereddened colors. 
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TABLE 7 
Catalog with Revised Teff 







SDSS 






IRFM'' 






KIC 








KIC ID 


Tcff 




0"raii 


Tcff 


0"lot 




Tcff 




[Fe/H] 




flag'= 




(K) 


(K) 


(K) 


(K) 


(K) 


(K) 


(K) 


(dex) 


(dex) 


(K) 




757076 


5137 


85 


55 


5150 


98 


94 


5174 


3.60 


-0.08 








757099 


5523 


97 


34 


5270 


110 


101 


5589 


3.82 


-0.21 








757137 


4822 


74 


42 


4536 


101 


99 


4879 


2.58 


-0.08 


49 





757218 


4728 


79 


17 


4489 


90 


75 


4555 


2.28 


-0.12 


67 





757231 


4909 


116 


64 


4974 


111 


89 


4825 


2.60 


-0.08 


24 






Note. — Effective temperatures presented here were computed at a fixed [Fe/H]= -0.2. 
^ Teff estimates based on / - A^^ using the original formula in CIO. 

^ Teff correction for giants. The sense is that this conection factor has been subtracted from the SDSS 7^.f|- estimate in the above table. 
^ Quality flag indicating stars with unusually discrepant SDSS Teif estimates (see text). 



• Use our griz-T^ff polynomials (Table |2|l or the origi- 
nal isochrone (Table[Tli to obtain temperature estimates. 
If complementary IRFM estimates are desired, use the 
CIO polynomials (for VtJKs). 

• Adjust hot-end temperatures above 6000 K using equa- 
tions 5-7. The polynomials in Table|2]are for the origi- 
nal SDSS temperature calibration (Table[TJ without the 
hot-end adjustment described in Section [33] 

• In Table |2] we adopted a metallicity [Fe/H]= -0.2 and 
a dispersion of 0.3 dex for error purposes. We also 
adopted a fractional error of 15% in the extinction. 

• The SDSS temperatures are inferred from the weighted 
average of the independent color estimates using the 
photometric errors discussed in Section |2l and the ran- 
dom uncertainties are the maximum of the formal ran- 
dom errors and the dispersion in those inferred from 
different griz colors. 

• If the metallicity is known independent of the KIC, the 
SDSS temperatures can be corrected using the values 
in Tableland if desired the IRFM temperatures can be 
corrected for metallicity by adopting star-by-star metal- 
licities in the CIO formulae. 

• Apply gravity corrections in Table |4] for giants with 
logg(KIC) < 3.5. 

• Outside the temperature range of the SDSS calibration, 
zero-point shifts of 223 K at the hot end and 150 K at 
the cool end should be applied to the KIC T^ff to avoid 
artificial discontinuities in the temperature scale at the 
edges of validity of the method. 

• In our revised Tgff table, we did not apply statistical cor- 
rections for binaries, but the current Table |6] could be 
employed to do so, and this should be included in pop- 
ulation studies. 

• We expect about 4% of the sample to have photome- 
try impacted by blends. Such stars could be identified 
as those having an excess dispersion from individual 
SDSS colors on the order of 130 K or more, and/or as 
those showing more than a 3cr deviation from the mean 
difference between the IRFM and SDSS temperatures. 



4.2. Main T^ff Catalog 

Our main result, the revised Teff for 1 61,9 77 stars in the 
long-cadence KIC, is presented in Table ITFl All of our re- 
vised Teff estimates in the catalog are based on the recalibrated 
griz magnitudes in the KIC (Section [2. 2l l. In addition to the 
^nz-based SDSS T^n, Table |7] contains (J - K,)-hased IRFM 
Teff using the original CIO relation, and KIC values along with 
logg and [Fe/H] in the KIC. The null values in the SDSS 
Teff column are those outside of the color range in the model 
(4043 K < Teff < 7509 K). Similarly, the CIO IRFM Teft are 
defined at 0.07 <{J- K,)o < 0.80. 

Statistical properties of our final temperature estimates are 
listed in Table [8] for dwarfs and for giants, separately. The 
relative KIC, IRFM, and SDSS temperatures for dwarfs and 
giants in the final catalog are compared in Figure [18] These 
comparisons include the adjustment to the hot end published 
SDSS scale described in Section [331 We did not correct the 
IRFM temperature estimates for gravity effects in the giants. 
The discrepancy between the two scales for the cool giants 
is consistent with being caused by this effect, as can be seen 
from the gravity sensitivity of (J - Ks) in Figure |7] 

Below we describe each column of Table |71 and provide a 
summary on how to correct Tgff for different logg, binarity 
(blending), and metalhcity. 

4.2. 1 . Error Estimates in Tc[[ 

For the SDSS and IRFM, we estimated total (ctot) and ran- 
dom ((Tian) errors for individual stars as follows. The random 
errors for the SDSS were taken from two approaches, tabu- 
lating whichever yields the larger value: a propagated error 
from the photometric precision and the one from measure- 
ments of Teff from individual color indices (g - r, g - i, and 
g - z). For the former, we repeated our procedures of solving 
for Teff with 0.01 mag photometric errors in gri and 0.03 mag 
errors in z: we added corresponding Teff errors from individ- 
ual determinations. The random errors for the IRFM were 
estimated from the 2MASS -reported photometric errors in J 
and Ks (combined in quadrature). 

In Table |2] we included systematic errors from ±15% er- 
ror in the foreground dust extinction and ±0.3 dex error in 
[Fe/H] from our fiducial case ([Fe/H]= -0.2) for both SDSS 
and IRFM measurements. The total error (titot) is a quadrature 
sum of both random and systematic error components. The to- 
tal errors are dominated by the extinction uncertainties, which 

' ' Only a portion of this table is shown here to demonstrate its form and 
content. A machine-readable version of the full table is available online. 
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TABLE 8 
Statistical Properties of T^^s 
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Note. — Statistical properties derived from the full long-cadence sample, after applying the hot 7[.fi conections. No metallicity and binary corrections were applied. 
^ Weighted mean difference {Tcii), weighted standard deviation (cr), and the expected dispersion propagated from random errors (cTprop). 
^ Median standard deviation of ^r/^-based temperature estimates from g - r, g - i, and g - z. 
^ Median dispersion expected from photometric err ors in griz. 



relate to both galactic position and distance. The quoted val- 
ues yield dispersions in temperature between YREC, IRFM, 
and spectroscopy consistent with the data. We present effec- 
tive temperatures defined at a fixed [Fe/H]= -0.2. If it is de- 
sired to correct for metallicities different from this fiducial 
[Fe/H], Teg corrections in Table[3]can be used. 

4.2.2. Corrections for different \ogg 

Our application of the isochrone assumes that all of the 
stars are main-sequence dwarfs. To correct for differences 
between the KIC and the model log^ values, we used log? 
sensitivities of the griz colors using Castell i & K urucz (2005 
ATLAS9 models, as described in Section 13.21 Table |4] lists 
the correction factors in T^s as a function of each color in- 
dex over Alog^ = 0.5-3.0 in a 0.5 dex increment. For a 
given color in each of these color indices, a difference be- 
tween the KIC and the model log g can be estimated (A log g = 
loggKic - loggYREc), and the corresponding ATeff values in 
Table |4] can be found in g - r, g - i, and g- z, respectively. 



The mean ATeff correction was then added to the dwarf-based 
Teff estimates. Our catalog (Table|7]) lists SDSS Tiff estimates 
already corrected using these \ogg corrections for those with 
logg(KlC) < 3.5 at 7;ft(SDSS) < 5300 K. If it is desired to 
recover the dwarf-based solution, correction terms (ATeff) in 
Table |2] should be subtracted from the listed reff(SDSS). 

4.2.3. Corrections for Binaries 

As described in Section [331 unresolved binaries and blend- 
ing can have an impact on the overall distribution of photo- 
metric Teff. If the population effect is of greater importance 
than individual Teff, correction factors in Table |6] should be 
added to the SDSS and IRFM (making them hotter) in Ta- 
ble |7] With l%-3% errors in griz photometry, it is difficult 
to distinguish between single stars with unresolved binaries 
and/or blended sources in the catalog. 

4.2.4. Quality Control Flag 
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Fig. 18. — Comparisons of Teff using the final SDSS Teff estimates. Comparisons are shown for the original KIC Teff for dwarfs (top left) and giants (top right), 
and for the (J — A^5)-based IRFM estimates for dwarfs (bottom left) and giants (bottom right). 



The last column in Table |7] shows a quality control flag. If 
the flag is set(flag= l),the SDSS Tgff values should be taken 
with care. The flag was set 

• if the standard deviation of individual T^<s from three 
color indices (g - r, g - i, and g - z) exceeds 130 K 
(A^= 1,402) 

• or if the difference between SDSS and IRFM measure- 
ments is greater than 3<t random errors (summed in 
quadrature) with respect to the mean trend (A^ = 4, 388). 
Only those at 4700 < Teff < 7000 K for dwarfs and 
4700 < Teff < 5400 K for giants were flagged this way, 
to avoid a biased ATeft distribution at the cool and hot 
temperature range (see Figure [THIl. 

• or if any of the griz measurements are not reported in 
the KIC (A^ = 257). 

In total, 5,798 stai's (about 4% of 154,931 stars with a valid 
SDSS Teff) were flagged this way. 

4.3. IRFM T^s from Tycho-2MASS System 

In addition to our main catalog in Table [T] we also present 
in Table|9]the IRFM T^n in Tycho Vj and 2MASS JHK, colors 
for 7,912 stars. These stars are a subset of the long-cadence 
KIC sample, which are bright enough to have VV magnitudes, 
and can be used as an independent check on our Teff scale (see 
lower left panel in Figure |6]l. The IRFM Teff values are pre- 
sented using Vj - J,Vt - H, Vt - Ks, and J - K^, with both 



random (CTian) and total (utot) errors. As in Table|7] random er- 
rors are propagated from photometric uncertainties, and total 
errors are a quadrature sum of random and systematic errors 
(15% error in reddening and 0.3 dex error in [Fe/H]). 

5. SUMMARY AND FUTURE DIRECTIONS 

The Kepler mission has a rich variety of applications, all 
of which are aided by better knowledge of the fundamental 
stellar properties. We have focused on the effective temper- 
ature scale, which is a well-posed problem with the existing 
photometry. However, in addition to the revised KIC temper- 
ature there are two significant independent results from our 
investigation. We have identified a modest color-dependent 
offset between the KIC and SDSS DR8 photometry, whose 
origin should be investigated. Applying the relevant correc- 
tions to the KIC photometry significantly improves the in- 
ternal consistency of temperature estimates. We have also 
verified that the independent temperature scales (Johnson- 
Cousins and SDSS) of An et al. and those from recent IRFM 
studies (Casagrande et al.) are in good agreement, permitting 
a cross-calibration of the latter to the SDSS filter system. Be- 
low we summarize our main results for the KIC, then turn to 
the major limitations of our main catalog, a brief discussion 
of the imphcations, and prospects for future improvements. 

5.1. Summary 

Our main result is a shift to higher effective temperatures 
than those included in the existing KIC. We have employed 
multiple diagnostic tools, including two distinct photometric 
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Note. — Effective temperatures presented here were computed at a fixed [Fe/H]= -0.2. 



scales and some high-resolution spectroscopy. In the case 
of cool (below 6000 K) dwarfs, the various methods for as- 
signing effective temperature have an encouraging degree of 
consistency. The Johnson-Cousins measurements of An et al. 
(l2007a ) are in good agreement with the independent IRFM 
temperatures from CIO in star clusters. In Table|5] for exam- 
ple, the V - Ic results agree within 15 K for all clusters if we 
adopt the Sandquist (2004) dataset for M67. The SDSS-based 
A09 s ystem is c onstructed to be on the same absolute scale 
as the lAn et al.l (,2007a) system, so a similar level of agree- 
ment is expected between the IRFM and the temperatures that 
we derive from the SDSS filters. A comparison of the IRFM 
and SDSS temperatures in the KIC confirms this pattern, with 
agreement to better than 100 K for the cool stars. Even this 
level of disagreement overestimates the underlying accord in 
the systems, because the IRFM (7 - Ks) diagnostic that was 
available to us in the KIC has systematic offsets relative to 
other IRFM thermometers even in the open clusters. When we 
correct for these offsets, the agreement for cool stars between 
the SDSS-based method of A09 and the IRFM (7 - K,) tem- 
peratures is very good, with average differences below 25 K 
and maximum differences below the 50 K level. Our cool 
dwarf temperatures are also within 50 K on average when 
compared to the spectroscopic results from Bll. The spectro- 
scopic sample of MZl 1 is cooler at the 88 K level, which we 
take as a measure of systematic uncertainties in the spectro- 
scopic scale (See Bruntt et al. 2010, for a further comparison 
of the spectroscopic and fundamental temperature scales). 

For hotter dwarfs the revised temperature estimates are 
higher than in the KIC, but the magnitude of the offset is 
not consistent between the two photometric scales and the 
spectroscopic data. Motivated by this offset, we adjusted the 
SDSS-based system of A09 to be cooler on average by 100 
K between 6000 K and 7000 K on the IRFM system. The 
consistency between photometric and spectroscopic scales de- 
grades for stars in this range. This could reflect defects in the 
fundamental temperature scale for hotter stars; the existing 
fundamental data for the IRFM include relatively few solar- 
abundance dwarfs above 6000 K. There could also be errors 
in photometric or spectroscopic temperature estimates from 
the onset of rapid rotation above 6300 K, or color anomalies 
from chemically peculiar hot stars. On the spectroscopic side, 
it would be valuable to compare the atmospheric temperatures 
inferred from Boltzman and Saha constraints to fundamental 
ones; as discussed in CIO, there can be significant systematic 
offsets between these scales for some systems. This issue de- 
serves future scrutiny and additional fundamental data would 
be very helpful. 

In the case of evolved stars we also found a hotter tem- 
perature scale than in the KIC. We had to employ theoretical 



estimates of gravity sensitivity, however, to temperature diag- 
nostics derived for dwarfs. An extension of the fundamental 
work to giants has been performed for other colors in the past, 
and it would be beneficial to test the theoretical predictions 
against actual radius data. 

5.2. Cautions and Caveats in Usage of the Catalog 

There are some significant drawbacks of the existing cata- 
log, and care is required in its proper application. Binary com- 
panions will modify the colors and temperatures of stars; we 
have provided tables for statistical corrections, but have not 
included this in the tabulated effective temperatures. Blend- 
ing can also impact colors, and there is clear evidence of some 
blended objects in our comparison of the KIC to SDSS DR8 
data with superior resolution. The major error source for the 
temperature estimates is the uncertainty in the extinction. We 
have adopted a global percentage value based on typical er- 
rors in extinction maps, but there could be larger local vari- 
ations. The color combinations available to us have limited 
diagnostic power for star-by-star extinction and binary cor- 
rections. For population studies, the stars in the long-cadence 
KIC sample were selected for a planet transit survey, and do 
not represent an unbiased set of the underlying population. 

The KIC abundance estimates have significant errors, 
largely because the filters with the greatest metallicity sen- 
sitivity were not available. As a result, we have adopted 
metallicity insensitive temperature diagnostics, but the tem- 
peratures should be corrected for individual metallicities if 
available. These effects are at the 100 K per dex level, and 
will therefore be smaller than the extinction uncertainties for 
most stars in the sample. The \ogg values for hot stars are not 
well-constrained in the KIC, but we have adopted KIC gravi- 
ties for cool stars. Our results would be affected at a modest 
level by changes in the derived gravities, and the appropriate 
corrections should be made if precise values are available. 

There are two open areas for further discussion as well: the 
appropriate temperature scale for the hot dwarfs and errors in 
the photometry. In the former case, we recommend adjust- 
ments above 6000 K to the SDSS scale. For the entire domain 
we also note inconsistencies between the J - Ks calibration 
and the other color-temperature relationships in the IRFM. 
Even after putting the fundamental photometric temperature 
scale on a common system, however, there is a difference be- 
tween it and the spectroscopic scale for stars above 6000 K. 
Until it is resolved we recommend inclusion of systematic 
temperature errors in this domain. The impact of the \ogg 
determinations on the extinction estimates for the hot stars 
should be investigated as well. The gravity diagnostics for the 
hot stars are not well measured, and asteroseismic gravities 
confirm this expected lack of precision. The KIC catalog in- 
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eluded this as an ingredient in the distance estimates, but it 
is difficult to reconstruct the weights and importance of this 
uncertainty after the fact. Star-by-star extinctions would be 
useful for this purpose. 

The origin of the differences between the SDSS (DR8) and 
KIC photometry should also be tracked down, and there may 
be spatially dependent or magnitude dependent terms. We 
also noted some cases with severe internal inconsistency in 
the photometric temperature diagnostics, and flagged those 
which we identified. We believe that unresolved blends are 
a promising candidate, but further work on this front is war- 
ranted. In a small fraction of cases these photometric issues 
can cause severe errors in the temperatures. Effective tem- 
peratures for stars where different colors return very different 
estimates should be treated with caution. 

Despite these reservations, we believe that the addition of 
temperatures more closely tied to the fundamental scale will 
significantly improve the reliability of inferences about the 
underlying stellar populations. 



5.3. Implications and Future Directions 

A shift to higher effective temperatures will have conse- 
quences for both planetary and stellar science. On the main 
sequence, hotter stars will be on average more massive and 
larger This would imply larger planet radii on average for 
such objects. The radii of evolved stars require more infor- 
mation (especially from surface gravity effects), and the con- 
sequences of the temperature scale shift for them are more 
difficult to predict from first principles. Stars of known as- 
teroseismic radius will be on average more luminous, which 
could partially explain dis crepancies in the mas s-radius rela- 
tionship for evolved stars jChaplin et al. II20TTI) . Asteroseis- 
mic parameters defined with scaling relationships will also be 
impacted. A more precise absolute effective temperature scale 
will also permit more stringent constraints on asteroseismic 
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pro perties from detailed modeling of the frequency spectrum 
(see iMetcalfe et afllMoh . 

However, the full potential will be realized as comple- 
mentary information becomes available on the Kepler sam- 
ple. Blue data (such as Johnson UB or: SDSS u) could be 
employed to infer more reliable photometric metallicities; 
Johnson-Cousins UBV(RI)c data would enable more reliable 
extinction estimates, binary discrimination, and broader ap- 
plication of the IRFM directly to stars in the sample. Pho- 
tometric systems naturally designed for F-type stars, such as 
Stromgren, would be useful for addressing the temperature 
and surface gravity scales in that regime. 

A more robust set of input data would provide an impor- 
tant control sample for the measured planet population; it will 
be challenging to obtain spectroscopic temperatures of both 
the planet candidates and the background stellar population. 
A better calibration of the fundamental temperature scale is 
possible once asteroseismic radii are combined with paral- 
laxes in the Kepler field, either via Kepler data or through 
the Gaia mission. The time domain data from the satellite are 
exquisite; a proper application of complementary tools from 
stellar astrophysics is now essential to fully realize their con- 
siderable scientific promise. 
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